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Abstract. The simplest modeling of planar quantum waveguides is the Dirichlet eigenproblem for the Laplace 
operator in unbounded open sets which are uniformly thin in one direction. Here we consider V-shaped guides. 
Their spectral properties depend essentially on a sole parameter, the opening of the V. The free energy band is a 
semi-infinite interval bounded from below. As soon as the V is not flat, there are bound states below the free energy 
band. There are a finite number of them, depending on the opening. This number tends to infinity as the opening 
tends to (sharply bent V). In this situation, the eigenfunctions concentrate and become self-similar. In contrast, 
when the opening gets large (almost flat V), the eigenfunctions spread and enjoy a different self-similar structure. 
We explain all these facts and illustrate them by numerical simulations. 



Introduction 

A quantum waveguide refer to nanoscale electronic device with a wire or thin surface shape. In the first case, 
one speaks of a quantum wire. The electronic density is low enough to allow a modeling of the system by a 
simple one-body Schrodinger operator with potential 

1 — > -Atp + Vip in M^. 

The structure of the device causes the potential to be very large outside and very small inside the device. As a 
relevant approximation, we can consider that the potential is zero in the device and infinite outside ; this can be 
described by a Dirichlet operator 

-0 I — > —Alp in J7 and t/j = on 

where U is the open set filled by the device. We refer to [3] where we can see, at least on the numerical 
simulations, the analogy between a problem with a confining potential and a Dirichlet condition. 

These kinds of device are intended to drive electronic fluxes. But their shape may capture some bound states, 
i.e. eigenpairs of the Dirichlet problem: 

—AtP = \ip in Q. and V = on (90. 

The topic of this paper is two-dimensional wire shaped structures, i.e. structures which coincide with strips of 
the form x (0, a) outside a ball of center and radius R lai^ge enough. These structures can be called planar 
waveguides. More specifically, there are the bent waveguides and the broken waveguides: Bent waveguides have 
a constant width around some central smooth curve, see Figure 1 , and the central curve of a broken waveguide 
is a broken line, see Figure 2. 

Due to the semi-infinite strips contained in such waveguide, the spectrum of the Laplacian — A with Dirichlet 
conditions is not discrete: It contains a semi-infinite interval of the form [p, +oo) which is the energy band 
where electronic transport can occur. The presence of discrete spectrum at lower energy levels is not obvious, 
but nevertheless, frequent. 

A remarkable result by Duclos and Exner [10] (and generalized in [6]) tells us that if the mid-line of a planar 
waveguide is smooth and straight outside a compact set, then there exists bound states as soon as the line is not 
straight everywhere. For broken guides, a similar result holds, [11, 2]: There exist bound states as soon as the 
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Figure 1. Curved guide Figure 2. Broken guide 

guide is not a straight strip. The question of the increasing number of such bound states for sharply bent broken 
waveguides has been answered in [5]. 

Before stating the main results of this paper, let us mention a few other works involving broken waveguides 
such as [15, 16]. It also turns out that our analysis leads to the investigation of triangles with a shaip angle 
as in [14] (see also [4]). The generalization to dimension 3 of the broken strip arises in [12] where a conical 
waveguide is studied, whereas a Born-Oppenheimer approach is in progress in [28]. 

In this paper, we revisit the results of [2, 5] and prove several other quantitative or qualitative properties 
of the eigenpairs of planar broken waveguides. Here are the contents of the present work: In section 1 we 
recall from the literature notions of unbounded self-adjoint operators, discrete and essential spectrum, Rayleigh 
quotients. After proving that the Rayleigh quotients ai^e increasing functions of the opening angle 9 of the 
guide (section 3), we adapt the technique of [6] to give a self-contained proof of the existence of bound states 
(section 4). We prove that the number of bound states is always finite, though depending on the opening angle 9 
(section 5), that this number tends to infinity like the inverse of the opening angle when ^ — )• (section 7). 
Concerning eigenvectors we prove that they satisfy an even symmetry property with respect to the symmetry 
axis of the guide (section 2). Their decay along the semi-infinite straight parts of the guide can be precisely 
evaluated (section 6) and is stronger and stronger when 9 decreases. We perform numerical computations 
by the finite element method, which clearly illustrate these decay properties. When the opening gets small, 
concentration and self-similarity appears, which can be explained by a semi-classical analysis: We give some 
overview of the asymptotic expansions established in our other paper [9]. We end this work by evaluations of 
the numerical convergence of the algorithms used for our finite element computations (section 9). 

Notation. The 1? norm on an open set U will be denoted by || • \\u- 

1. Unbounded self- adjoint operators 

In this section we recall from the literature some definitions and fundamental facts on unbounded self- 
adjoint operators and their spectrum. We quote the standard book of Reed and Simon [30, Chapter VIII] and, 
for the readers who can read french, the book of Levy-Bruhl [23] (see in particular Chapter 10 on unbounded 
operators). 

Let H he. & separable Hilbert space with scalar- product (•,•)//. We will consider operators A defined on a 
dense subspace Dom(A) of H called the domain of A. The adjoint A* of A is the operator defined as follows: 

(i) The domain Dom{A*) is the space of the elements uof H such that the form 

Dom(A) 3 V ^ {u,Av)h 
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can be extended to a continuous form on H. 
(ii) For any u G Dom(yl*), A*u is the unique element of H provided by the Riesz theorem such that 

(1) Vi; G Dom(A), {A*u,v)h = {u,Av)h 

Definition 1.1. The operator A with domain Dom(A) is said self-adjoint if A = A*, which means that 

(2) Dom(yl*) = Dom(A) and Vu, -u G Dom(A), {Au,v)h = {u, Av)h- 

1.1. Operators in variational form. The operators that we will consider can be defined by variational formu- 
lation. Let us introduce the general framework first. Let be given two separable Hilbert spaces H and V with 
continuous embedding of V into H and such that V is dense in H. Let b be an hermitian sesquilinear form on 
V 

b:V xV 3 {u,v)^ b{u,v) G C 

which is assumed to be continuous and coercive: This means that there exist three real numbers c, C and A 
such that 

(3) Vu G y, c\\u\\v <b{u,u) + A{u,u)h < C\\u\\y . 
Let A be the operator defined from V into its dual V' by the natural expression 

Vf G V, {Au,v)h = b{u,v). 

In other words, for all u £ V, Au is the linear form v i— )• b{u, v). Note that the operator A + A Id is associated 
in the same way to the sesquilinear form 

b + A(-, ■)h --V xV 3 {u,v) ^ b{u, v) + A(m, f )h G C 

which is strongly elliptic by (3). As a consequence of the Riesz theorem (or the more general Lax-Milgram 
theorem) there holds 

(4) ^ + AId is an isomorphism from V onto V' . 

This situation provides many examples of (unbounded) self-adjoint operators. The following lemma is re- 
lated to the Friedrichs's lemma cf. [23, Section 10.7]. 

Lemma 1.2. Let A be the operator associated with an hermitian sesquilinear form b coercive on V. Let A be 
the restriction of the operator A on the domain 

Dom(A) = {ueV : AueH]. 

Then A is self-adjoint. 

Proof. The operator A is symmetric because b is hermitian. Moreover we check immediately that 

yu,v G Dom(A), {Au,v)h = {u,Av)h = b{u,v). 

In particular, we deduce Dom(A) C Dom(A*). Let us prove that Dom(A*) C Dom(yl). Since the domain of 
A and A* are, unchanged by the addition of A Id, and since 

(^ + Aid)* = A* + Aid 

we can consider ^4 + A Id instead of A, or, in other words, using (4), assume that A is bijective. 

Then we deduce that A is an isomorphism from Dom{A) onto H: Indeed, A is injective because A is 
injective; If / G there exists u £ V such that Au = f, and u G Dom{A) by the very definition of Dom(A). 

Let w belong to Dom(A*). This means, cf. (2), that w £ H and A*w = f £ H. Let us prove that w belongs 
to Dom(A). 
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Since A is bijective there exists u G Dom(A) such that Au = f. We have 

Vf G y, b{u,v) = {f,v)H 

and therefore 

Vt; G Dom(yl), {u,Av)h = {'w,Av)h- 
Hence, as A is bijective, for all g ^ H, {u, g)H = {w, g)H- Finally u = w, which ends the proof. □ 

Exemple 1.3. Let Q. be an open set in and do\M a part of its boundary. On 

y = {^eH^(J7): V = OonaDir^^} 

we consider the bilinear form 

b{il),i)') = [ V^(x) • VV''(x) dx. 

The operator A is equal to — A and it is self-adjoint on H = L'^{Q,) with domain 

Dom(^) = {V' G y : AipGh'^in) and a„V = on 50 \ aoir!^}. 

1 .2. Discrete and essential spectrum. Let A be an unbounded self-adjoint operator on H with domain Dom (A) . 
We recall the following characterizations of its spectrum a{A), its essential spectrum aess{A) and its discrete 
spectrum ad\s{A): 

• Specti'um: A G cr{A) if and only if {A — Aid) is not invertible from Dom (A) onto H, 

• Essential spectrum: A G cress(^) if and only if {A — X Id) is not Fredholm^ from Dom(^) into H (see 
[30, Chapter VI] and [23, Chapter 3]), 

• Discrete spectrum: ad\s{A) := a{A) \ aessiA). 

We list now several fundamental properties of essential and discrete spectrum. 

Lemma 1.4 (Weyl criterion). We have A G aess{A) if and only if there exists a sequence (un) G Dom(A) such 
that 1 1 tin II// = 1, (un) has no subsequence converging in H and {A — AId)n„ — )• in H. 

From this lemma, one can deduce (see [23, Proposition 2.21 and Proposition 3.11]): 

Lemma 1.5. The discrete spectrum is formed by isolated eigenvalues of finite multiplicity. 

Lemma 1.6. The essential spectrum is stable under any perturbation which is compact from Dom (A) into H. 

Exemple 1.7. Let A be the self-adjoint operator on H associated with an hermitian sesquilinear form b coercive 
on V, cf. Lemma 1.2. Let us assume that V is compactly embedded in H. Then the spectrum of A is discrete 
and formed by a non-decreasing sequence Uk of eigenvalues which tend to +oo as A; — )• +oo (see [23, Chapter 
13]). Let (ufc)fc>i be an associated orthonormal basis of eigenvectors: 

Avk = UkVk, Mk>l. 

Then we have the following identities 

(5) WuGH, \\u\\h = '^\{u,Vk)jjf, 

k>l 

(6) G y, b{u,u) = y^^Uk\{u,Vk)fff, 

k>l 

(7) VnGDom(A), WAuWl = Y,^l\{u,Vk) 

k>l 

^We recall that an operator is said to be Fredholm if its kernel is finite dimensional, its range is closed and with finite codimension. 
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Exemple 1.8. Let us define A^"^ as the Dirichlet problem for tfie Laplace operator —A on an open set C M", 
cf. Example 1.3 with Soir^ = 90.. 

( i) If Q is bounded, A^"^ has a purely discrete spectrum, which is an increasing sequence of positive numbers. 

( ii) Let us assume that there is a compact set K such that 

0, \ K = \^ Oj (disjoint union) 

j finite 

where Qj is isometrically affine to a half-tube = (0, +oo) x toj, with loj bounded open set in M"~^. Let fij 
be the first eigenvalue of the Dirichlet problem for the Laplace operator — A on coj. Then, we have: 

<Tess(AQ"') = Uj[;Uj,+oo) = [min/ij, +oo). 

j 

The proof can be organized in two main steps. Firstly, for each j we construct Weyl sequences supported in 
Sj associated with any A > fij, which proves that cress(AQ"') C [miiij Hj, +oo). Secondly we apply Lemma 
1.6 with A = B — C where A = A^" and B = A^"' + W, where M/^ is a non negative and smooth potential 
which is compactly supported and such that W > miiij fij on K. On one hand, since C is compact, we get that 
(Jess(^) = o'essiB). On the other hand, we notice that: 

/ llVV'f dx+ [ VF|V'pdx> [ llV^f dx + min/ij / IV'P dx 
Jn Jn Jq i Jk 

and, using the Poincare inequality with respect to the transversal variable in each strip : 

/ llV^f dx > V / llV^f dx > V / dx > mmfij [ \^|;\^ 

Jn ~^ JT.j ~^ JT.j 3 Jn\K 



We infer that: 



J J^j ■ JT.j Jn\i 

[ llVV'f dx+ [ dx > min^j [ l^/^p dx. 

Jn Jn i Jn 



The min-max principle provides that inf (j{B) > min^ fij so that inf (Tess(-B) > minj fij, and finally we get 

info"ess(^) ^ min^j. 

j 

For an example of this technique, we refer for instance to [6, Section 3.1]. Let us notice that the Persson's 
theorem provides a direct proof (see [29] and [13, Appendix B]). 

The same formula holds even if the boundary conditions on dQ n K are mixed Dirichlet-Neumann. 

1.3. Rayleigh quotients. We recall now the definition of the Rayleigh quotients of a seLf-adjoint operator A 
(see [23, Proposition 6.17 and 13.1]). 

Definition 1.9. The Rayleigh quotients associated to the self-adjoint operator ^ on i7 of domain Dom(A) are 
defined for all positive natural number j by 

Aj = mi sup — — . 

Mi,...,UjeDom(A) we[iii,.--,"j] V^T^'^IH 
independent 

Here [ui, . . . , Uj] denotes the subspace generated by the j independent vectors ui, . . . , Uj. 

The following statement gives the relation between Rayleigh quotients and eigenvalues. 

Tlieorem 1.10. Let Abe a self-adjoint operator of domain Dom(A). We assume that A is semi-bounded from 
below, i.e., there exists A G M such that 

V'uGDom(A), {Au,u)h + ^u,u)h 
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We set 7 = minc7ess(^)- Then the Ray lei gh quotients Xj of A form a non-decreasing sequence and there holds 

(i) If Xj < 7, it is an eigenvalue of A, 

(ii) If Xj > 7, then Xj = 7, 

(Hi) The j-th eigenvalue < 7 of A (if exists) coincides with Xj. 

Lemma 1.11. Let A be the self-adjoint operator on H associated with an hermitian sesquilinear form b coer- 
cive on V, cf. Lemma 1.2. Then the Rayleigh quotients of A are equal to 



u-i,...,Uj£V uelui,...,Uj] \U,U)h 
independent 

Corollary 1.12. Let A and A be the self-adjoint operators on H and H associated with the hermitian sesquilin- 
ear forms b and b coercive on V and V, respectively. We assume that 

H C H, V CV, b{u, u) > b{u, u) Vu G V. 

Let Xj and Xj be the Rayleigh quotients associated to A and A, respectively. Then 

Vj > 1, Xj > Xj . 

Exemple 1.13 (Conforming Galerkin projection). This consists in choosing a finite dimensional subspace V 
of V, which also defines a subspace of H, and b = b. The Rayleigh quotients Xj are the eigenvalues of the 
discrete operator A, and they are larger than the Rayleigh quotients Xj of A. 

Exemple 1.14 (Monotonicity of Dirichlet eigenvalues). Let $7 be an open subset of M" and 17 C il. The 
extension by from to realizes a natural embedding of HQ(il) into HQ(il), and of L^(il) into L^(0). 
Then, with obvious notations: 

Vj > 1, Xji^f) > A,(A8'^). 

Remark 1.15 (Non-monotonicity of Neumann eigenvalues). The Neumann problem consists in taking H^{Q,) 
as variational space. The argument above does not work because there is no canonical embedding of H^{Q,) 
into H^{Q). Moreover, the monotonicity with respect to the domain is wrong: 

( i) Let us choose Q bounded and connected, and Q a subset of with two connected components. Then 

Ai(A^-) = A2(A^-) = and AilA^ ) = 0, X^iA^^ > 0. 

(ii) If we take two embedded intervals for Q, and Q,, then explicit calculations show that Aj(A^^") > Aj(A[^^"). 

Another nice and non trivial counter-example can be found with the de Gennes operator appearing in the 
superconductivity theory, cf. [8]. 

2. The broken guide 

Let us denote the Cartesian coordinates in by x = (xi,X2). The open sets U that we consider are 
unbounded plane V-shaped sets. The question of interest is the presence and the properties of bound states for 
the Laplace operator A = df -\- 82 with Dirichlet boundary conditions in such il. We can assume without loss 
of generality that our set Q is normalized so that 

• it has its non-convex corner at the origin = (0, 0), 

• it is symmetric with respect to the xi axis, 

• its thickness is equal to vr. 
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The sole remaining parameter is the opening of the V: We denote by ^ G (0) f ) the half-opening and by Jig the 
associated broken guide, see Figure 3. We have 

(8) = ((xi,X2) G : xi tan6' < |x2| < (xi + ) tan6' I . 

I V sm uy ) 

Finally, like in Example 1.8, we specify the positive Dirichlet Laplacian by the notation A[^^^ This operator is 
an unbounded self-adjoint operator with domain 




Figure 3. The broken guide VIq (here 9 = |). 



The boundary of is not smooth, it is polygonal. The presence of the non-convex corner with vertex at 
the origin is the reason for the domain Dom(AQ|^'^) to be distinct from n Yi\{Q,0). Nevertheless this domain 
can be precisely characterized as follows. Let us introduce polar coordinates (p, ip) centered at the origin, with 
9? = coinciding with the upper part X2 = xi ianO of the boundary of Q.q. Let x be a smooth radial cutoff 
function with support in the region xi tan0 < \x2\ and x = 1 in a neighborhood of the origin. We introduce 
the explicit singular function 

(9) V'sing(2;i,X2) = x(p)/9"/''sin^, with u = 2{Ti-e). 
Then there holds, cf. the classical references [22, 18]: 

(10) Dom(AH;) = h2 n ^l{^e) © [V'sing] 
where [V'sing] denotes the space generated by V'sing- 



2.1. Essential spectrum. 

Proposition 2.1. For any 6 € (0, ^) the essential spectrum of the operator A[^^'' coincides with [1, +c«). 

Proof. This proposition is a consequence of Example 1.8 (ii): Outside a compact set, Q.q is the union of two 
strips isometric to (0, +oo) x (0, tt). Since the first eigenvalue of —dy on Hq(0, tt) is 1, the proposition is 
proved. □ 
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2.2. Symmetry. In our quest of bounded states (A, '(/') of Aq|^^ i.e. 

I — = in 0,0, 



(11) 



with A < 1, V / 0, 



ip = on d^le 
we can reduce to the half-guide Q,'^ defined as (see Figure 4) 

~ {(^1' ^2) € 0,e '■ X2 > 0}, with the Dirichlet part of its boundary doir^Q^ = dilg n dQ^, 

Q 



as we ai^e going to explain now. Let us introduce Aj^l'^ as the positive Laplacian with mixed Dirichlet-Neumann 



conditions on 0^ with domain, cf. Example 1.3: 

Dom{A^P = {tl; e H^{nj) : At/' G ^^(0+), = on SoirJ^^ and ^zV" = on xa = O}. 

Then (7ess(A^|'^) coincides with aessi^Qg)- Concerning the discrete spectrum we have: 



(-ii^'O) 




-^1 



Figure 4. The waveguide Qg and the half-guide (here 9 = |). 
Proposition 2.2. For any 6 G (0, f ), adis(A^p coincides with adU^^p- 

Proof. The proof reUes on the fact that A^J^'' commutes with the symmetry S : (xi, 2:2) ^ {xi, — X2). 

(i) If (A, Tp) is an eigenpair of Aj^+ , the even extension of ux to defines an eigenfunction of A^"' associated 
with the same eigenvalue A. Therefore we have the inclusion 

adis(Af;ix)Cadis(AB:). 

(ii) Conversely, let (A, ^p) be an eigenpair of A^^'' with A < 1. Splitting ip into its odd and even parts ip"^^ and 
^even ^jjj^ rcspcct to X2, wc obtain: 

^J; = V^°dd ^ ^even^ A^^"'''' = A^°^^ and A^^"^"" = At/;^^^". 

We note that tp"^^ satisfies the Dirichlet condition on the line X2 = 0, and ip^^^" the Neumann condition on the 
same line. Let us check that xp°^^ = 0. If it is not the case, this would mean that A is an eigenvalue for the 
Dirichlet Laplacian on the half- waveguide J]^. But, by monotonicity of the Dirichlet spectrum with the respect 
to the domain, cf. Example 1.14, we obtain that the spectrum on Q'^ is higher than the spectrum on the infinite 
strip which coincides with Q,^ when xi > 0. This latter spectrum is equal to [1, +00). Therefore the Dirichlet 
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Laplacian on the half- waveguide $7^ cannot have an eigenvalue below 1. Thus, we have necessarily: tp"^^ = 
and i/j = i/j^^^" which is an eigenfunction of A'^|'^ associated with A. □ 

We take advantage of Proposition 2.2 for further proofs and for numerical simulations. 

3. MONOTONICITY OF RAYLEIGH QUOTIENTS WITH RESPECT TO THE OPENING 

From now on we consider the Rayleigh quotients associated with the operator A^+, which we write in the 
form, cf. Lemma 1.11 



(12) Xj{e) = inf sup 

tpi,...,ipj independent in Hi(rj+), '/'6[i/'i v.'/'j] 

lt!l,...,ll!j=0 on doir^'^ 

Those Xj{0) which are < 1 are all the eigenvalues of A[^|J sitting below its essential spectrum. 

Proposition 3.1. For any integer j > 1, the function i— )• Xj{9) defined in (12) is non-decreasing fi-om (0, ^) 
into M+. 

Proof. We cannot use directly Corollary 1.12 because of the part of the boundary where Neumann conditions 
are prescribed. Instead we introduce the open set VIq isometric to Q.'^ , see Figure 5, 

= G (~:[^'+°°) ^ (0,7r) : y < xtan6l + 7r if x G (^--^,o)| . 

(0,^) 



y = 7r 



-> (-vr,0)'^ y = 

Figure 5. The reference half-guide Q, := ^t^/^. 

The part d^-ir^e of the boundary carrying the Dirichlet condition is the union of its horizontal parts. The 
numbers Xj{6) can be equivalently defined by the Rayleigh quotients (12) on fig. 
Let us now perform the change of variable: 

X = X tan 6, y = y, 

so that the new integration domain Q := i7^/4 is independent of 6. The bilinear- gradient form b on is 
transformed into the anisotropic form be on the fixed set O: 

(13) be{i>,^') = ftan^9{d,ijd,i;') + i^y^P ^y^p') dxdy, 



with associated form domain 

(14) V := {^p e R\n) : iP = on dow^} 

independent of 9. 

The function 9 i— )• tan^ 9 being increasing on (0, ^), we have 

y-ip ^ V, 9 bg{'ip, V) non-decreasing on (0, §). 
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We conclude thanks to Corollary 1.12. □ 

Remark 3.2. Using the perturbation theory, we also know that, for all j > 1, the function 6 i— )• Xj{9) is 
continuous with respect to G (^if)- Moreover, Ai being simple (as the first eigenvalue of a Laplace- 
Dirichlet problem), it is analytic because we are in the situation of an analytic family of type (B) (see [21, p. 
387 and 395]). In fact the numerical simulations lead to think that all the eigenvalues below 1 are simple and 
thus analytic. 



We recall that the lower bound of the essential spectrum of the operator Aj^l'^ is 1. Its first Rayleigh quotient 



4. Existence of discrete spectrum 

, Mix 

is given by, cf. (12), 

I|VVII^+ 

(15) Ai(^)= inf 

In this section, we prove the following proposition: 

Proposition 4.1. For any 6 G (0, the first Rayleigh quotient Xi{9) is < 1. 

This statement implies that Xi{9) is an eigenvalue of A^!^ (application of Theorem 1.10), hence of the 

Dirichlet Laplacian A^^ on the broken guide (Proposition 2.2). 

This statement was first established in [2]. Here we present a distinct, more synthetic proof, using the method 
of [6, p. 104-105]. 

Proof. For convenience, it is easier to work in the reference set Cl = Clj^/^ introduced in the previous section, 
with the bilinear form bg (13) and the form domain V (14). We are going to work with the shifted bilinear form 



be{ip,ip') = b0{'4),'il)') - li)ij'dxdy 

Jn 

which is associated with the quadratic form: 

QeW = be{^P,iJ). 

Then the first Rayleigh quotient of Qg is equal to Xi{0) — 1. To prove our statement, this is enough to construct 
a function tf; ^ V such that: 

QeW < 0. 

This will be done by the construction of 

(1) A sequence ipn such that Qe{ilJn) — )■ as n — oo, 

(2) An element cf) oiV such that bg {ipn , ip) is nonzero and independent of n. 

The desired function will then be obtained as a suitable combination ipn + e(f). Let us give details now. 

Step 1. In order to do that, we consider the Weyl sequence defined as follows. Let x be a smooth cutoff 
function equal to 1 for x < and for x > 1. We let, forn G N \ {0}: 

Xnix) = X (^) and ^n{x,y) = Xn{x)sm.y. 
Using the support of Xn, we find that Qeii/^n) is equal to 

/O pX+TT poo pIT 

/ (cos^ y - sin^ y) dy dx + / / ( tan^ 6'(x^)^ sin^ y + x^(cos^ 2/ - sin^ y)) dy dx. 
'TtJo Jo Jo 
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Then, elementary computations provide: 

I (cos^ y — s\v? y)dy = and / / (cos^ y — sin^ y) dy dx = 0. 

Jo J-TT Jo 



Moreover, we have: 



X iu)\ du ] . 



roc rn / rl 

/ / tan2e(x;)2sin2ydydx< / 
Jo Jo \Jo 

Hence we have proved that Qe{ipn) tends to as n — )• oo: 

(16) QeiiPn) < ^ with Ke = (^j^ \x'{u)\^ du^ vrtan^ 9 . 

Step 2. We introduce a smooth cutoff function r] of x supported in (— vr, 0). We consider a function / of 
y € [0,vr] to be determined later and satisfying /(O) = 0. We define (p{x,y) = r]{x)f{y). For e > to be 
chosen small enough, we introduce: 

V'n,e(a;, y) = iJn{x, y) + e(t){x, y). 

We have: 

Let us compute be{Tpn, <P)- We can write, thanks to considerations of support: 

/O r-X+TT r-0 pX+tt 

/ ri{x){cosyf'{y) - sin yf{y)) dydx = / r]{x){cosyf{y))' dydx. 
-TT Jo J-TT Jo 

Using /(O) = 0, this leads to: 

/■o 

heiil^nA) = I cos(x + 7r)/(x + vr) dx. 

J —TT 

We choose f{y) = 7]{y — vr) cos(y — vr) and we find: 



This implies, using (16): 



fO 

bei'^n, 0) = — / '>f{x) cos^(x) dx = — F < 0. 



Qe{'^n,e) <^-2Te + De\ 



where D = Qg{4>) is independent of e and n. There exists e > such that: 

-2Fe + De^ < -Fe. 
The angle 9 being fixed, we can take N large enough so that 

2N - 2^' 

from which we deduce that Qe{'^N,£) < —eV/2 < 0, which ends the proof of Proposition 4.1. □ 
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5. Finite number of eigenvalues below the essential spectrum 



For a self-adjoint operator A and a chosen real number A we denote by M{A, A) the maximal index j such 
that the j-th Rayleigh quotient of A is < A. By extension of notation, if the operator A is defined by a coercive 
hermitian form 6 on a form domain V, and if Q denotes the associated quadratic form Q{u) = b{u, u), we also 
denote by J\f{Q, A) the number M{A, A). This is coherent with the fact that in this case the Rayleigh quotients 
can be defined directly by Q, cf. Lemma 1.11: 



This section is devoted to the proof of the following proposition: 

Proposition 5.1. For any 6 G (0, |), 7V(A[^j^'', 1) is finite. 

Thus in any case has a nonzero finite number of eigenvalues under its essential spectrum. 

Proof. For the proof of Proposition 5. 1 we use a similar method as [27, Theorem 2. 1]. 

Like for the proof of Proposition 4.1, it is easier to work in the reference set (l introduced in section 3, with 
the bilinear form hg (13) and the form domain F (14). The opening 9 being fixed, we drop the index 9 in the 
notation of quadratic forms and write simply as Q the quadratic form associated with hs: 



We recall that the form domain V is the subspace of -i/; G H^{^) which satisfy the Dirichlet condition on do\rU.. 
We want to prove that 

N{Q, 1) is finite. 
We consider a partition of unity (xo; Xi) such that 

Xo{xf +Xi{xf = 1 

with Xo{x) = 1 for X < 1 and Xo{x) = for x > 2. For i? > and i e {0, 1}, we inti^oduce: 

xiA^) = xi{R~^x). 

Thanks to the IMS formula (see for instance [7]), we can split the quadratic form as: 



inf sup 

U\,...,Uj(iV u^[ui,...,u 

independent 



Q{u) 

ij] {u, u)h ' 




(17) 



IIX0,R^ll| - \\Xl,Rll^% 



We can write 



X0,Ri^)\' + \Xl,R{^)\' = R 



^Wr{x) with 



Wr{x) = \x'o{R~'x)\^ + \xi{R-'x)\^ . 



Then 



(18) 




Let us introduce the subsets of $7: 



Oo,R = {{x,y) €n 



X < 



2R] 



and Oi^R = {(x, y) e ^ : x > R] 
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and the associated form domains 

Vo = ^(pe H^{Oo,r) : = on dow^^ D dOo,B^ and on {2R} x (0, 7r)| 

Vi = RI{Oi,r). 
We define the two quadratic forms Qo,r and Qi^r by 

(19) Qi,R{^)= I tan" e\d,<P\^ + \dycl)\^ - R-^WRixM^dxdy for ^P G Ve, £ = 0,1. 
As a consequence of (17) and (18) we find 

(20) QW = QoMXO,Ri^) + Ql,R{Xl,Ri^) V^GF. 
Let us prove 

Lemma 5.2. We have: 

AfiQ,!) <Af{Qo,R,l)+Af{Qi,R,l). 
Proof. We recall the formula for the j-th Rayleigh quotient of Q: 

Xj = mf sup --— 2" . 

ECV i,&E M\~ 
dim E=j " 

The idea is now to give a lower bound for Xj. Let us introduce: 

V ^ VoxVi 

As {xo,r,Xi,r) is a partition of the unity, J is injective. In particulai\ we notice that J : V ^ JiY) is 
bijective so that we have: 

Xj = mf sup --— 2" 

F^J(V) ^,6^-i(i.) IIV'II^ 
dim F=j 



J: 



sup — n TiTp n 

FCJ{V) ^^j-i^F) \\XQ,R'^\\~ + Xl.flV' 5 
dimF=j ^ ^ ' ^ ^ 



Qo,r{XO,R'^) + Ql,RiX l,Ri') 
llxo.ijV'll- + \\Xl,R^ 

= mt sup ,, , II , • 

FCjjV) (^o,^,)6F IIV-oll^^^ + IIV'lllo,^ 
dim F=j 

As J^{V) C Vb X Vi, we deduce by an application of Corollary 1.12: 
(21) > ml sup II , ||„ II , 112 =: Vj, 



dimF 



Let Ag^R be the self-adjoint operator with domain Dom(^^ /j) associated with the coercive bilinear form corre- 
sponding to the quadratic form Qe,R^ on Ve. We see that Uj in (21) is the j-th Rayleigh quotient of the diagonal 
self-adjoint operator Ar 

i^O^ A^ ) with domain Dom{Ao R) x Dom{Ai^R) . 
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The Rayleigh quotients of A^ ji are associated with the quadratic form Qe^R for £ = 0,1. Thus vj is the j-th 
element of the ordered set 

{Afc(Qo,i?), A: > 1} U {Xk{Qi,R), k > 1}. 
Lemma 5.2 follows. □ 

The operator Aq^r is elliptic on a bounded open set, hence has a compact resolvent. Therefore we get: 
Lemma 5.3. For all R > 0, M{Qo^r, 1) is finite. 

To achieve the proof of Proposition 5.1, it remains to establish the following lemma: 
Lemma 5.4. There exists Rq > such that, for R > Rq, M^Qi^r, 1) is finite. 
Proof. For all G Vi, we write: 

(j) = Uo(l) + Uicj), 

where 

(22) IlQ4>{x,y) = ^{x)shiy with <I>(x) = / (j){x , y) sin y dy 

Jo 

is the projection on the first eigenvector of —dy on Hj!|(0, vr), and Hi = Id — Ho. We have, for all e > 0: 

QiM^) = QiA^Q^) + QiA^i^) - 2 / R^^WR{x)nQ(l) Ui^dx dy 
(23) 

> Qi,ij(no0) + Qi,i?(ni,/.) - / R~^WR{x)\Uocl)\^dxdy-e [ R-^WR{x)\Uicj)\^ dx dy 

Since the second eigenvalue of —dy on Hq(0, vr) is 4, we have: 

[ \dyUicl)\^dxdy>mim^^. 
Denoting by M the maximum of Wr (which is independent of R), and using (19) we deduce 

QiM'^i^) > i^-MR-^m.m^^^. 

Combining this with (23) where we take e = 1, and with the definition (22) of IIo, we find 

Qi,r{<^) > QRi'^) + (4 - 2MR-^)\\n,m^^^, 

where 

g^($) = / tan^ eid^^l"^ + l^l"^ - R''^Wr{x)\^\'^ dx 
Jr 

POO 

> / tan'^e\d^^\'^ + 1^1"^ - R-^Ml[R2R]\^\'^dx. 
Jr 

We choose R = y/M so that (4 - 2MR-'^) = 2, and then 

(24) >g/?W + 2||ni0||2,^^, 
where now 

POO 

(25) qR{^) = / tan^ 9\d,^\^ + (1 - 1[r.2/?])|^|' dx. 

JR 
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Let CLR denote the ID operator associated with the quadratic form q^. From (24)-(25), we deduce that the j-th 
Rayleigh quotient of Ai^r admits as lower bound the j-th Rayleigh quotient of the diagonal operator: 

CLR 

2Id 



so that we find: 



N{Ql,B.^)<^f{qB.^)■ 



Finally, the eigenvalues < 1 of clr can be computed explicitly and this is an elementary exercise to deduce that 
M{qB, 1) is finite. □ 



This concludes the proof of Proposition 5.1. 



□ 



6. Decay of eigenvectors at infinity - Computations for large angles 

6.1. Decay at infinity. In order to study theoretical properties of eigenvectors of the operator A[^JJ correspond- 
ing to eigenvalues below 1, we use the equivalent configuration on Q.q introduced in section 3, see also Figure 
6. The eigenvalues < 1 of A^^*^ ai^e the same as those of A-"* (with Dirichlet conditions on the horizontal parts 

of the boundary of Q.q) and the eigenvectors are isometric. The main result of this section is a quasi-optimal 
decay in the straight part (0, +oo) x (0, tt) of the set as a; — ;> oo. 



(0,^) 



(- 



vr 



X 



tan 6 



0) 



y = 7r 



y = 



Figure 6. The half-guide (here 6 = f ). 



Proposition 6.1. Let 9 G (0, ^). Let ip be an eigenvector o/ A-"^ associated with an eigenvalue A < 1. Then 
for all e > the following integral is finite: 

(26) 



^0 



Proof. We give here an elementary proof based on the representation of Tp as solution of the Dirichlet problem 
in the half-strip S := R+ x (0, vr) 



(27) 



-Alp = Xip in S, 

tpix, 0) = 0, V(5 , vr) = Vx > 0, 
^(0,y) = <7 VyG(0,^) 
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where g is the trace of ip on the segment F, see Figure 6. Since if) belongs to H^(S), its trace on 9S belongs 
to H^/^. Because of the zero trace on the lines y = and y = vr, we find that g belongs to the smaller space^ 
Hqq^(F), which is the interpolation space of index \ between Hj(F) and L2(F), cf. [24]. 

We expand ijj on the eigenvector basis of the operator —dy, self-adjoint on n Hq(0, tt). Its normalized 
eigenvectors are 

Vk{y) = sin ky with eigenvalue z^^. = k'^. 

We expand g in this basis: 

9{y) = y23k'"kiy), where 5^=/ 9iy)vk{y)dy. 
k>i -^0 

Interpolating between (5) and (6), we find 

k>l 

We can easily solve (27) by separation of variables. We find 



(28) ^Pii, y) = Yl e-^v^ gk Vk{y). 

k>l 

The estimate of ^|; in (26) is then trivial. Let us prove now the estimate of dxip- We use that X]fc>i ^ 9k finite 
so that we can write: 

= -Y^y^k^e-'''^gkVk{y). 

k>l 

We have: 

gx(v^-e)5_^(~^ ~) = -5]\/fc^e*(^-^^) e-'^gkVk{y) 

k>l 

leading to the estimate: 

poo poo 

/ / |e^(v^-)5,^(x,y)fdxdy = / {k'-X)e''i^'-^')e-'^'gldi 

Jo Jo fc>l-^0 



oo 



<(ykgl)snp fce-2^^-^^e-2-dx, 
^ fc>i Jo 



where 7 = 7(A) > is a constant, uniform with respect to A; > 1. Using the change of variables x 1— > X, 
we can see that the integrals 



are uniformly bounded as k — 00, which ends the proof of the estimate of dxip in (26). The estimate of dyi/j is 
similar-. □ 



^ The space HQQ^(r) is the subspace of H^'''^(r) spanned by the functions v such that 

y{'^-y) 
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Remark 6.2. Under the conditions of Proposition 6.1, we have also a sharp global L°° estimate 

||e*v^VllL^(n,) <oo. 

1/2 

To prove this we use again the representation (28) and the fact that the trace g is more regular than Hqq (F). In 
fact, as a consequence of the characterization (10) of the domain of the operator ^q^, the trace g belongs to 
HQ(r) and, therefore, the sum Y2k>i 9k finite by (6). Then, we have: 



k>l 



Taking absolute values, we get: 



-El 



1/2 



k>l 



1/2 



fc>l fc>l 

Remark 6.3. The estimate (26) can also be proved with a general method due to Agmon (see for instance [1]). 

6.2. Computations for large angles. When — ^, Xii9) tends to 1, see [2]. The representation (28) shows 
that in such a situation, the behavior of the associated eigenvector t/j is dominated by its first term, proportional 
to 

e~ 



-xVi^ sin({/). 




Artificial boundary 



Figure 7. The model half-guide := 



In order to compute such an eigenpair by a finite element method, we have to be careful and take lai^ge 
enough domains — we simply put Dirichlet conditions on an artificial boundary far enough from the corners 
of the guide. Our computations are performed in the model half-guide Q := ^^^^^ for the scaled operator 



(29) 



^9 



equivalent to — A in through the variable change 

sin 9 and 



u 



X2 



\f2cos9. 



We use a Galerkin discretization by finite elements in a truncated subset of Q. with Dirichlet condition on the 
artificial boundary, see Figure 7. According to Corollary 1.12, cf. also Examples 1.13 and 1.14, the eigenvalues 
)^^^{9) of the discretized problem are larger than the Rayleigh quotients \j{9) of ^g. When the discretization 

gets finer and the computational domain, lai^ger, X^^if)) tends to \j{9) for j = 1,2,... 
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e = 0.6889 * tt/2 XT\^) = 0.98844 



Figure 8. Computations for moderately large angles. Plots of the first eigenvectors in the 
physical domain (rotated by |). Numerical value of the corresponding eigenvalue Xi{6). 

For the values of 9 considered in this section (9 > j), the numerical evidence is that the discrete spectrum of 
^0 has only one element Ai(^). The numerical effect of this is the convergence to 1 of all other computational 
eigenvalues XY^{9) for j = 2,3, . . . 

The computations represented in Figure 8 are performed with the artificial boundary set at the abscissa 
u = 5tt\/2. The plots are mapped back to the corresponding physical domain by a postprocessing of the 
numerical results. 

The computations in Figure 9 are performed with the artificial boundary set at the abscissa u = 107r\/2. The 
plots use the computational domain because the corresponding physical domains would be too much elongated 
to be represented. 

7. Accumulation of eigenpairs for small angles 

When 9 tends to 0, there is more and more room for eigenvectors between the two corners of the guide. 
For any rectangular box B contained in like in Figure 10, by the monotonicity of Dirichlet eigenvalues 
(Example 1.13), we know that for any j 

X,{9) < X,{B) 
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6" = 0.7022 *7r/2 6* = 0.8538 * 7r/2 6^ = 0.9702 * 7r/2 

* (61) = 0.9903037 Af* (61) = 0.9994215 A^p* (6* ) = 0.9999998 



Figure 9. Computations for very large angles with the mesh M4L (see Figure 17). Plots in 
the computational domain i7. 




Figure 10. The waveguide and a Dirichlet box B inside. 

where Xj{6) are the Rayleigh quotients of A^'^ and Xj{B) the Dirichlet eigenvalues on B. We choose B in the 
form of a rectangle bounded by the vertical lines xi = — avr and xi = 0, and the horizontal lines X2 = ±/37r. 
Thus Q and /3 satisfy, c/ (8) 

1 vr 
aG(0, ^— ), /37r = (-asin6' + 1)- 



sm a cos d 

The eigenvalues of the Dirichlet problem in B are 

^2 £2 

We look for the eigenvalues Xj (B) less than 1. Therefore a has to be chosen > 1. Thus /3 < (1 — sin 6) / cos 6 < 
1 for any 9. As a consequence k = 1 and the eigenvalues Aj(;B) less than 1 are necessaiily of the form 

cos^ 7^ 
+ — 



4(1 -a sin 6^)2 

We optimize B: The minimum of Xj{B) is obtained for a such that 

sin 9 cos^ 9 if' 
2(1 - a sin 61)3 ~ ~ 
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Figure 11. Computations for 6 = 0.0226 * 7r/2 ~ 2° with the mesh M64S (see Figure 18). 
Numerical values of the 10 eigenvalues \j{6) < 1. Plots of the associated eigenvectors in the 
physical domain. 



Since we are interested in the behavior as ^ — 0, we take without asymptotic loss 
which provides 



cos^ d 



4V(l_4l/3j2/3 sin2/3^)2 

As a consequence, as soon as the quantity Z := 4^/3j2/3 gjj^^/s q less that the first root of the equation 
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i.e. for Z < 0.4679, we have Xj{B) < 1 and, hence, Xj{9) < 1. This imphes that the maximal number J such 
that \j{B) < 1 is greater than 

0.4679^/^ • 0.5 sin"^ 9 ~ 0.1601 sin~^ 9 . 

Therefore the number of eigenvalues of A^'^ less than 1 tends to infinity (at least) like as ^ tends to 0. 

We present in Figure 1 1 the computations of the all eigenpairs of A^'^ for the angle 9 = 0.0226 * 7:/2. We 
find 10 eigenvalues < 1. Note that for this angle 9, the numerical value of 0.1601 sin~^ 9 is 4.5103. 

8. Asymptotic behavior of eigenpairs for small angles 

We present in Figures 12 and 13 computations of the first eigenvector for smaller and smaller values of the 
angle 9. We notice that the eigenvectors look similar, with the appearance of a short scale in the horizontal 
variable. 




6* = 0.1482 *7r/2 6^ = 0.1032 * 7r/2 6* = 0.0701 * 7r/2 

A^P*(^) = 0.56209 Af*(^) = 0.48754 X'l^\9) = 0.42763 



Figure 12. Computations for small angles. Plots in the computational domain 

We perform in [9] asymptotic expansions of the first eigenpairs at any order as 6* — )■ 0, using techniques of 
semi-classical analysis (see for instance [20, 19] and also [26]). We briefly describe now some of the results 
proved there. 

Let us recall that the eigenvalues Aj (9) < 1 of our operator Aq^^ coincide with those of the scaled operator 

:= -2 sm'^9 dl - 2 cos'^9 dl 

in the model half-guide Q.. The construction and validation of asymptotic expansions for the eigenpairs of J^q 
rely on a Bom-Oppenheimer approximation. Roughly, this consists of four steps: 

(1) Definition of an operator -Sfg^^ in the horizontal variable u by replacing the operator in the vertical 
variable '■= —2cos^9 in each slice u = const, by its first eigenvalue A(n) 

(2) Semi-classical analysis of the eigenpairs of as ^ — ^ 0. 

(3) Determination of a change of variables {u, v) ^ {u, t) on in order to exhibit a tensor product 
structure. Here appears the role of the limit operator ^„ as n 0. Its first eigenvector is t; i— 
cos(f /2). 

(4) Construction of expansions of eigenvectors in the new variables. 
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(9 = 0.0416 *7r/2 6* = 0.0270 * 7r/2 = 0.0112 * 7r/2 

A^P*(0) = 0.37085 Af*(^) = 0.33845 X'l^\e) = 0.29766 



Figure 13. Computations for very small angles. Plots in the computational domain 0. 



Step 1: Definition of 

Let T{u) be the intersection of Q with the vertical line of abscissa u. 

• For u G (— 7r\/2, 0), X(n) = (0, n + 7r\/2). The operator has Neumann condition at and Dirichlet 
at n + Thus 



A(u) = 2cos26'' 



7r2 



4(n + 7r^/2)2 

• For u G (0, +oo), X(n) = {u,u + nV^). The operator ^ has Dirichlet conditions at both ends. Thus 

A(n) = cos^6'. 

The Bom-Oppenheimer operator is 
Step 2: Semi-classical analysis of 

The operator can be viewed as a ID Schrodinger operator with potential A. The potential has a well with 
bottom at n = 0. The well is not smooth and has a triangular shape, see Figure 14. 

The behavior of the eigenpairs of -Sfg^^ as ^ — is governed by the Taylor expansion of the potential A at 
the well bottom u = 0, i.e. by the tangent potential V defined by 

{1 u 
^, if tt < 0, 
4 2ttV2 
1, ifu>0. 

The con^esponding model problem is the problem of the behavior as /i — )• of the eigenpairs of the operator 

" |l, ifn>0. 

The potential barrier on (0, +oo) produces a Dirichlet condition at u = for the leading terms of the asymp- 
totics: We are led to the Airy-type eigenvalue problem 

(30) - h^dliph - uiPh = Ehi^h on (-oo, 0), with V(0) = 0. 
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Function A 
Function V 



Figure 14. The Born-Oppenheimer potential A and its left tangent at u = (here cos^0 is 
set to 1). 

The change of variables X = uh~'^/^ transforms the above equation into the reverse Airy equation 

-(9|-^ - X-^ = on (0, +oo), with l'(O) = 0, (where = h~^^^Eh) 

whose solutions can be easily exhibited: As we will see, the eigenvalues fi are the zeros of the reverse Airy 
function A(X) := A\(—X), where Ai is the standard Airy function. The zeros of A form an increasing sequence 
of positive numbers, which we denote by z/\{j), j > 1. 
Since -A"{X) - XA{X) = we have for any G M 

-A"{X) - {X - E)A{X) = EA{X) i.e. - f^' {X + E) - XA{X + E) = EA{X + E) . 

If E = zj\{j), then A{X + E) vanishes at X = 0, hence 

{z^{j),A[X + za(j))) is an eigenpair. 

Conversely, all eigenpairs are of this form. 

We deduce that the eigenvalues of problem (30) are E^ = h?^^z/^{j) and the associated eigenvectors are 

V^(u) =A(n/i-2/3 + ^^(j-)). 

Coming back to our operator ^g^^, we prove in[9] that its eigenvalues have asymptotic expansions of the 
form 

where 7^^^ are some real coefficients. The eigenvectors have expansions in powers of 01/3 

using the scale 

ii/i~2/3 for n < and n/i"^ for u > 0. 
Step 3: Tensorial structure of J^g. 

On the left part of Q, i.e. its triangular part Tri in the half-plane u < 0, we perform a change of variables to 
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transform Tri into a square: 

(32) Tri 9 (u, v) ^ {u, t) G (-7r\/2, 0) x (0, TrV2), with t = 

u + Try 2 

On the right part of Q., i.e. its strip part contained in the half -plane u > 0, we perform a change of variables to 
transform it into a horizontal strip: 

(33) (n, H> (n, r) G (0,+oo) X (0, 7rA/2), with t = v - u. 

This allows to work with the Taylor expansions of the transformed operator around u = 0, on the left and on 




Figure 15. The model waveguide $7 and the change of variables. 

the right. The operator = — 2c)| appears on the left with Dirichlet condition at t = 1 and Neumann at 
i = 0. Its first eigenvector is cos |, associated with the eigenvalue \. The operator ^_|. = —25^ appears on 
the right with Dirichlet condition at t = 0, 1. 

Step 4: Asymptotics of the eigenpairs of ^q. 

The projection on the eigenvector 1 1— )• cos | appears naturally at the first step of the expansion (this projection is 
sometimes called Feshbach or Grushin projection) and lets appear ■^^'^ . A complete asymptotics for eigenpairs 
can be constructed, and in a further step, validated. As a result \j{9) has an expansion similar to \^'^{9), with 
different coefficients forn > 3: 

where 7j „ are some real coefficients. In Figure 16 we display the functions {6 * 2/7r)^/'^ i— \j{0) for j = 
1 , 2 (computed by finite elements) and their linear approximation as ^ — >• 0, corresponding to the two-term 

approximation in (34), i.e. {9 * Ijixfl'^ ^ i + §3f^' 

In the reference domain 0, the eigenvectors of .ifg have a multiscale expansion in powers of 9^1'^. On the 
left side, there are two scales, the ultra-short scale {u9~^,v) and the short scale {u9~'^^^ ,v), on the right side 
there is only the ultra-short scale. The asymptotics is dominated by its first term, which is 

^ 2a(1))cos^. 



Q2/3 
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1st computed eigenvalue 
2d computed eigenvalue 
Essential spectrum 

Linear approximation of 1st eigenvalue 
Linear approximation of 2d eigenvalue 



0.1 0.2 0.3 0.4 0.5 0.6 



0.8 0.9 



Figure 16. The eigenvalues Ai and A2 as functions of {6 * 2/7r)^/^. 



This structure explains the results seen in Figures 12-13. 



9. Convergence of the finite element method 

We now present some aspects of the computation process that led to the numerical results shown in the 
previous sections. 

As described in the section 6.2, the eigenvalue problem writes ^gtp = Xtp in the domain 17 with homoge- 
neous Dirichlet conditions on its boundary, except on the horizontal segment where Neumann conditions are 
set, see identity (29) and Figure 7. The associate bilinear form is 



2 sin^ 6'(a„^ duij') + 2 cos^ e{dy7p dyip') du dv 



defined on the corresponding form domain 

V = {^p eU^in) : V = on ^Dirf^}- 
The eigenvalue problem writes in variational form: find non-zero ^ G y and A such that 

By Galerkin projection on a finite dimension subspace Vfds of V, this problem can be rewritten as the general- 
ized eigenvalue problem: find the eigenpairs (A, w) such that Sw = XMw, where S and M are the stiffness 
and mass matrices associated with a basis (^1, . . . , "^n) of Vf^s- 

and M=({^j,^k) 



V /i<i,fe<Af ' ^ y''^)/ i<j,k<N 

The computation process consists of two main steps: first, using a finite element method that leads to the two 
matrices S and M, and second, using an algorithm to compute the eigenpairs. In the following two sections, we 
focus on the algorithm for the computation of the eigenpairs and then on the influence of the choice of meshes 
and polynomial degrees in the finite element method. 
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All the computations have been done in double precision arithmetic, on a iMac computer (4 GB memory, 
3.06 GHz Intel Core 2 Duo processor). 

9.1. Computation of the eigenpairs. We have first written a program using the Fortran 77 finite element 
library "Melina" [25], which also provides a routine to compute eigenpairs of real symmetric eigenvalue prob- 
lems. This routine uses an algorithm based on subspaces iterations. We call Mel this method. 

We have also written another program using the C++ version of the previous library, called "Melina++". 
The computation of the eigenpairs have been made using the well-known package ARPACK++ [17]. We call 
Arp this method. We thus were able to reproduce the results obtained with the Mel method, which is a cross 
validation of both methods. 

Here, we compare the computation of the eigenpairs with the two methods. For this purpose, in both cases, 
we fix the parameters governing the finite element part: interpolation degree 6 at Gauss-Lobatto points, quad- 
rature rule of degree 13. 

We have selected two test configurations: 

• "large angle" configuration: 6 = 0.9702 * 7r/2, 4 eigenvalues computed, mesh M4L with 656 triangles 
(see Figure 17), leading to matrices of size 12325x12325; 

• "small angle" configuration: 9 = 0.0226 * 7r/2, 12 eigenvalues computed, mesh M64S with 6144 
triangles (see Figure 18), leading to matrices of size 111265x111265. 

Although the internal algorithms of the two methods are different, the parameters governing the computation 
of the eigenvalues are the same: a tolerance e that controls the end of the iteration process, and the dimension 
A'^sub of the subspaces involved in the subspace iteration — Nsuh is at least equal to Aval + 1> where A^^,al is 
the number of desired eigenvalues. We let A'sub vary between 10 and 70 and ran the programs for e = 10^", 
n = 4, 5, 6, 7, while recording the number of iterations and the CPU time needed. The CPU time of the finite 
element part, i.e. the computation of the matrices S and M, does not depend on these pai^ameters. It appeai^s 
on the graphs as CPUef. 

The results are gathered on Figure 19 for the "large angle" configuration and on Figure 20 for the "small 
angle" configuration. We can observe that the number of iterations decreases as N^uh increases. A horizontal 
line at the beginning of the first two graphs indicates a failed computation (no convergence after the chosen 
maximal number of iterations). For the values tested, the parameter e does not play a significant role on the 
number of iterations, nor really on the CPU time (although it does on the residual). Since the algorithms used 
in the two methods are not the same, the number of iterations cannot be compai^ed directly. They ai^e mainly a 
good indicator of the computation process behavior. 

It is also worth to notice that the memory requirements increase as the subspace dimension Ajub increases. 
This parameter is difficult to handle since it is problem dependent. These graphs suggest that there is no critical 
value: it can be chosen large enough to ensure computation to succeed, but not too large, mainly because of 
memory considerations. 

The CPU graphs of the Arp method seem a bit chaotic: both algorithms need an initial vector which is chosen 
as a random vector by default. From our experience, we can say that this method is more sensible to this initial 
vector than the Mel method. Finally, these graphs show clearly that the Arp method is more efficient than the 
Mel method. 

Remark 9. 1. We are interested in the smallest eigenvalues of the problem Sw = XMw. Both methods Arp and 
Mel provide an option to choose in which end of the spectrum the wanted eigenvalues are to be searched. Let 
us mention that these algorithms perform better at computing the largest eigenvalues in general (this is due to 
the fact that they are ultimately based on the power method). In our case, this is typical and computation nearly 
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Figure 17. Mesh M4L for large angle, 656 triangles. 




Figure 18. Mesh M64S for small angle, 6144 triangles. 
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30 40 50 

Subspace dimension 



70 




30 40 50 

Subspace dimension 



70 



Figure 19. Comparison of the behavior of the methods for the large angle 9 = 0.9702 * tt/2. 
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Figure 20. Comparison of the behavior of the methods for the small angle 6 = 0.0226 * 7r/2. 



always fails if the algorithms are asked to compute the smallest eigenvalues of this problem. For example, with 
ARPACK, we can observe the following behavior: 

• in the "large angle" configuration, for e = 10""* and A'sub = 17 or 45, the computation fails ; 

• in the "small angle" configuration, for e = 10~^ and Nsuh = 27 or 43, the computation fails ; if we 
change the mesh to a coarser one (M16S) with 384 triangles, the computation fails for iVsub = 27 and 
succeeds for Nsub = 43 with a CPU time of 236 sec. 
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Thus, the coiTect strategy is to compute the largest eigenvalues of the eigenvalue problem Mw = uSw and 
retrieve the wanted eigenvalues A by inversion {X = 1/u and the eigenvectors are the same). All the computa- 
tions presented here have been carried on using this strategy. To conclude this remark, let us mention that the 
computation that took 236 sec. with the wrong strategy, takes 0.47 sec. with the right one. 

9.2. Influence of meshes and polynomial degrees in the finite element method. We now choose the small 
angle 9 = 0.0226 * 7r/2 and fix the parameters governing the computation of the eigenvalues: e = 10~^, 
-^sub = 25, A^vai = 10, since there are 10 eigenvalues < 1 as shown on Figure 11. 

We have built five nested meshes of the same computational domain Q., called M4S, M8S (see Figure 21), 
M16S, M32S (see Figure 22) and M64S (see Figure 18). The number n in the name (MnS) is the number 
of segments of the subdivision of the horizontal (and vertical) boundary of Q. It indicates the chai^acteristic 
diameter h of the triangles, which is halved from a mesh to the next one in the list. Thus, the number of 
triangles is multiplied by 4 from a mesh to the next one. 

We have computed the 10 eigenvalues for e = 10^^ using the finest mesh M64S and considered the first and 
last eigenvalues obtained as reference values. We denote them by A^^^*" and A^q^ 

For each mesh, we let the interpolation degree k vary from 1 to 6, while recording the differences | Ai — A^^^| 
and |Aio — Aj^q"! obtained for each degree. The results are gathered on Figure 23. Each point on the graphs 
correspond to the result of a computation. Points corresponding to the same mesh have been linked together 
by a line to make the graphs more readable ; the corresponding degree is written below. The graphs show the 
convergence of the first and last computed eigenvalues: 

• with respect to the interpolation degree for a given mesh ; 

• with respect to the mesh, for a given degree. 

The number of degrees of freedom (d.o.f.) of the problem, which is the dimension of the matrices, is roughly 
proportional to (kn)^, which explains the choice of the abscissa log|Q(A^)/2. The precision attained is about 
one order of magnitude better for the first eigenvalue Ai than for the last one Aiq. For the first eigenvalue and 
the coarser meshes, we observe a kind of super convergence for the small degrees, then a linear convergence. 
The convergence tends to be linear as the mesh becomes finer. For the last eigenvalue Aio, the convergence is 
mainly linear. 
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Figure 21. Meshes M4S (24 triangles) and M8S (96 triangles). 
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Figure 22. Meshes M16S (384 triangles) and M32S (1536 triangles). 

e = 0.0226 * II/2 - k^= 0.327833 = 0.0226 * it/2 - 1^^ = 0.981 955 




log.uj(Number of d.o.f.)/2 log^jj(Number of d.o.f.)/2 



Figure 23. Convergence towards Ai ~ 0.32783 and Aio ~ 0.98195. 

This phenomenon of super convergence is due to a sort of locking which takes place in coarser meshes: such 
meshes with low degree k are not able to capture the fine scale structure of the eigenvectors (see in Figure 13 
representations of the first eigenvector in the computational domain for small angles). 

The average rate of convergence with respect to is twice the rate of convergence with respect to h, in 
accordance with a well-known convergence result in finite elements [31]. Typically the rate is 1 in /i and 2 in 
k^^. These somewhat low rates are due to the singularity at the reentrant corner of Q,. This singularity comes 
from the Laplace singularity (9), and for 9 = 0.0226 * tt/2, the singularity exponent ^ is ~ 0.5057. 

Nevertheless, all our computations are accurate enough to display clearly the asymptotic behavior of the 
eigenpairs in the small angle and large angle Umits, see Figure 16. 
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